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A  three-dimensional  nonhydrostatic  numerical  model  is  used  to  study  the  generation  of  internal  waves 
by  the  barotropic  tidal  flow  over  a  steep  two-dimensional  ridge  in  an  ocean  with  strong  upper-ocean 
stratification.  The  process  is  examined  by  varying  topographic  width,  amplitude  of  the  barotropic  tide, 
and  stratification  at  three  ridge  heights.  The  results  show  that  a  large  amount  of  energy  is  converted 
from  the  barotropic  tide  to  the  baroclinic  wave  when  the  slope  parameter,  defined  as  the  ratio  of  the 
maximum  ridge  slope  to  the  maximum  wave  slope,  is  greater  than  1.  The  energy  flux  of  internal  waves 
can  be  normalized  by  the  vertical  integral  of  the  buoyancy  frequency  over  the  ridge  depths  and  the 
kinetic  energy  of  the  barotropic  tides  in  the  water  column.  A  relationship  between  the  normalized 
energy  flux  and  the  slope  parameter  is  derived.  The  normalized  energy  flux  reaches  a  constant  value 
independent  of  the  slope  parameter  when  the  slope  parameter  is  greater  than  1.5.  It  is  inferred  that 
internal  wave  generation  is  most  efficient  at  the  presence  of  strong  upper-ocean  stratification  over  a 
steep,  tall  ridge.  In  the  Luzon  Strait,  the  strength  of  the  shallow  thermocline  and  the  location  of  the 
Kuroshio  front  could  affect  generation  of  internal  solitary  waves  in  the  northern  South  China  Sea. 
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1.  Introduction 

Large-amplitude  internal  waves  are  often  observed  in  satellite 
remote  sensing  images  in  areas  where  the  barotropic  tides  are 
forced  to  flow  over  a  two-dimensional  ocean  ridge  or  the 
continental  slope  (Global  Ocean  Associates,  2004).  Recent  field 
observations  have  described  the  structure  of  internal  waves  in  the 
northern  South  China  Sea,  (e.g.,  Ramp  et  al.,  2004;  Yang  et  al., 
2004),  Hawaiian  Ridge  (e.g.,  Nash  et  al.,  2006),  Bay  of  Biscay  (e.g. 
Gerkema  et  al.,  2004),  and  Massachusetts  Bay  (e.g.  Scotti  et  al., 
2007),  to  name  a  few.  In  the  Luzon  Strait,  internal  waves  produced 
at  the  steep  ridges  evolve  into  large-amplitude  internal  solitary 
waves,  propagating  hundreds  of  kilometers  westward  into  the 
northern  South  China  Sea  (Zhao  and  Alford,  2006;  Chao  et  al., 
2007).  The  flow  velocity  and  thermocline  depression  could  reach 
2  m/s  and  150  m,  respectively  (Ramp  et  al.,  2004;  Yang  et  al., 
2004).  In  the  deep  ocean,  the  large  vertical  motion  associated  with 
internal  waves  is  of  concern  to  maritime  activities.  When 
propagating  shoreward,  breaking  internal  waves  dissipate  and 
lead  to  turbulent  mixing  (e.g.,  Mourn  et  al.,  2003),  affecting 
plankton  distribution,  pollution,  and  sediment  transport  (Butman 
et  al.  2006). 

The  baroclinic  energy  flux  F  associated  with  internal  wave 
generation  over  a  ridge  is  likely  affected  by  the  topographic  width  L, 
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topographic  height  h0,  total  water  depth  H,  amplitude  of  tidal 
current  H0,  forcing  frequency  co,  buoyancy  frequency  N,  and  the 
Coriolis  parameter  /.  The  nondimensional  parameters  associated 
with  these  quantities  are  reviewed  by  Garrett  and  Kunze  (2007). 
An  important  parameter  is  the  slope  parameter  y,  the  ratio  of  the 
topographic  slope  to  the  slope  of  the  internal  wave  beam 
[(co2-/2)/(N2-co2)]1/2.  A  ridge  is  classified  as  supercritical  if  y  >  1 
and  subcritical  if  y  <  1 .  Other  relevant  parameters  are  the  relative 
topographic  height  3  =  h0  /H,  the  local  Froude  number  Fr  =  U0 /( Nh0 ), 
and  the  tidal  excursion  parameter  RL  =  Uo /(coL). 

Linear  theory  of  internal  wave  generation  is  reviewed  by 
Garrett  and  Kunze  (2007).  The  solution  in  an  ocean  of  infinite 
depth  in  the  limit  y  1  is  given  by  Bell  (1975).  The  infinite-depth 
solution  establishes  a  scaling  factor  for  the  baroclinic  energy  flux 

F«p0U2/12(N2_w2)1/2(1_/2/w2)1/2j  (1) 

where  p0  is  a  reference  density.  The  h l  dependence  for  an  isolated 
Gaussian  bump  has  been  confirmed  by  Balmforth  et  al.  (2002). 
Eq.  (1)  has  been  used  by  nearly  all  theoretical  studies  on  internal 
wave  generation.  In  an  ocean  with  depth-varying  stratification, 
the  linear  theory  of  Llewellyn  Smith  and  Young  (2002)  suggests 
that  the  buoyancy  frequency  at  the  ridge  crest  should  be  used  for 
N  in  Eq.  (1),  and  the  wave  slope  is  determined  by  a  depth- 
averaged  N  using  the  WKB  theory. 

The  energy  flux  in  the  solution  of  Bell  (1975)  is  independent  of 
the  ridge  slope  (Llewellyn  Smith  and  Young,  2002).  The 
dependence  of  F  on  y  has  been  studied  analytically  by  several 
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investigators.  For  an  isolated  Gaussian  bump  in  an  infinitely  deep 
ocean,  Balmforth  et  al.  (2002)  found  that  F  is  14%  above  Bell’s 
(1975)  prediction  as  y  approaches  1.  The  prediction  of  Bell  (1975) 
could  be  significantly  reduced  when  the  ocean  is  of  finite  depth 
(Llewellyn  Smith  and  Young,  2002;  Khatiwala,  2003).  Of  more 
interest  is  the  range  y  >  1 .  In  the  extreme  case  of  a  knife-edge 
ridge  (7  ->  00)  with  small  h0,  the  energy  flux  is  exactly  twice  that 
for  a  bell-shaped  ridge  with  small  slope  (St  Laurent  et  al.,  2003; 
Llewellyn  Smith  and  Young,  2003).  Furthermore,  the  energy  flux 
increases  monotonically  as  the  relative  height  <5  increases  and 
reaches  two  times  Bell’s  (1975)  prediction  at  <5  =  0.92  (St  Laurent 
et  al.,  2003;  Llewellyn  Smith  and  Young,  2003).  Petrelis,  et  al. 
(2006)  showed  that  Bell’s  (1975)  solution  is  a  good  prediction  for 
subcritical  ridges  while  the  knife-edge  solution  provides  a  rough 
estimate  for  supercritical  ridges. 

At  the  critical  slope,  the  solution  is  singular  in  the  linear 
theory.  Wave  beams  appear  when  y  approaches  1  (Balmforth 
et  al.,  2002;  Khatiwala,  2003;  Petrelis  et  al.,  2006)  or  when  a  knife- 
edge  ridge  is  used  (St  Laurent  et  al.,  2003;  Llewellyn  Smith  and 
Young,  2003).  In  the  nonlinear  regime,  Khatiwala  (2003)  found 
that  energy  flux  of  the  internal  wave  may  saturate  for  certain 
topographic  shapes  over  a  supercritical  ridge.  Transient  hydraulic 
jumps  over  a  supercritical  ridge  are  demonstrated  in  the 
nonhydrostatic  numerical  simulation  of  Legg  and  Klymak 
(2008).  This  overturning  process  may  saturate  the  energy  flux  of 
large-amplitude  internal  waves  over  a  supercritical  ridge,  making 
the  energy  flux  less  dependent  on  the  slope  parameter. 

The  effect  of  RL  >  1  is  known  to  generate  internal  waves  of 
higher  tidal  harmonics  both  analytically  (Bell,  1975)  and  numeri¬ 
cally  (Legg  and  Huijts,  2006).  In  the  Luzon  Strait,  typical 
amplitude  of  the  barotropic  tidal  current  is  0.1  m/s,  and  the 
excursion  parameter  RL  is  much  less  than  one.  Thus,  energy  of 
internal  waves  in  the  Luzon  Strait  occurs  mostly  at  tidal 
frequencies.  The  Froude  number  is  small  over  the  tall  ridges 
reaching  the  thermocline  in  the  Luzon  Strait.  Nonlinearity  in  the 
Luzon  Strait  is  mainly  due  to  large  topographic  slope  and  the  large 
ratio  of  the  ridge  height  to  water  depth  <5.  In  this  study, 
nonhydrostatic  simulation  is  performed  to  investigate  the 
mechanisms  of  internal  wave  generation  over  a  tall,  steep  ridge 
in  an  ocean  with  strong  stratification  in  the  upper  ocean.  The 
internal  wave  energy  flux  is  calculated  to  find  its  dependence  on  y 
at  different  ridge  heights. 

The  paper  is  organized  as  follows.  Section  2  describes  the 
model  setup.  The  model  results  are  presented  in  Section  3  for 
three  sets  of  experiments  at  different  ridge  heights.  Effects  of  the 
topographic  width,  barotropic  tidal  amplitude,  and  stratification 
on  the  internal  wave  energy  flux  are  studied  in  each  set  of 
experiments.  A  parameterization  scheme  for  energy  flux  is 
derived  and  compared  with  order-of-magnitude  estimation  from 
the  real  ocean  in  Section  4,  followed  by  discussion  and  conclusion. 


2.  Model  description  and  experiments 

The  numerical  experiments  in  this  paper  extend  the  study  of 
Shaw  et  al.  (2009)  in  the  parameter  space.  The  three-dimensional 
nonhydrostatic  numerical  model  of  Shaw  and  Chao  (2006)  is  used. 
The  model  solves  the  three-dimensional  momentum,  continuity, 
and  density  equations  with  Boussinesq  and  rigid-lid  approxima¬ 
tions.  The  governing  equations  are 

Du  ,  1  0  ,  .  ^2  d2u 

-+2flk’xn  =  --Vp-^k+4VHu+vF  (2) 

(3) 


°£=KHVlP  +  K§.  (4) 

where  u  is  the  three-dimensional  velocity  vector  (u,  v ,  w)  in  the 
(x,  y,  z)  coordinates,  k  is  a  unit  vector  pointing  upward  from  the 
North  Pole,  k  is  the  local  upward  unit  vector,  p  is  the  deviation  of 
density  from  a  reference  density  po(1028  kg  m-3),  Q  is  the  Earth’s 
rotation  rate,  and  p  is  the  pressure.  The  kinematic  eddy  viscosity 
and  diffusivity  are  constant,  so  that  AH=  1 06  m2  s“ \  I<H=  1 05  m2  s“ \ 
v=\ x  10-4  m2  s-1,  and  k=0A  x  10-4  m2  s-1.  Eq.  (2)  includes 
both  the  horizontal  and  vertical  components  of  the  Coriolis 
acceleration.  The  finite  difference  form  of  the  governing  equations 
and  the  solution  procedure  are  described  in  Shaw  and  Chao 
(2006). 

The  model  domain  is  from  x=  -  200  to  200  km  in  the  east- west 
direction  and  from  z=  —H=  -800  m  at  the  bottom  to  0  m  at  the 
surface.  The  latitude  of  the  Luzon  Strait  (20°N)  is  used  in 
calculating  the  Coriolis  parameter.  The  horizontal  and  vertical 
grid  sizes  are  200  and  10  m.  In  the  y  direction,  two  grid  cells  are 
used.  With  a  periodic  boundary  condition,  no  y-variation  is 
introduced  if  the  initial  values  are  independent  of  y.  In  effect,  the 
three-dimensional  model  performs  a  two-dimensional  simulation 
in  this  study.  At  x=0  is  a  bell-shaped  (Witch-of-Agnesi)  ridge 

h(x)  =  /im/[l+(x/L)2]  (5) 

where  hm  is  the  maximum  height  and  L  is  the  half-width  of  the 
ridge.  The  maximum  ridge  slope,  located  at  x=31/2  L,  is  (33/2/8  )/im/L 
In  this  study,  hm  is  fixed  at  600  m,  and  L  is  varied  to  change  the 
maximum  slope.  Without  affecting  the  maximum  slope,  the  ridge  is 
displaced  vertically  to  alter  the  ridge  height  h0,  which  is  measured 
from  the  bottom  (z=  -800  m)  to  the  ridge  crest. 

Flow  over  the  ridge  is  forced  by  the  oscillating  barotropic  M2 
tide.  The  velocity  in  the  x-direction  is  given  by 

U  =  U0sin(cot),  (6) 

where  co  =  1/1 2.42  h-1  =  1.41  x  10“4s_1  is  the  tidal  frequency. 
The  barotropic  tidal  current  is  0  at  t=0  and  becomes  positive  in 
the  first  half  tidal  cycle.  In  this  paper,  positive  tidal  velocity  is 
called  ebb  tide  as  in  the  Luzon  Strait. 

The  perturbation  density  p'  is 

p'(x,z,t)  =  p-pa(z).  (7) 

The  ambient  stratification  pa(z)  is  given  by 

Pa(z)=~^  [l+tan/l/^)].  (8> 

where  A p  is  the  scale  of  density  perturbation.  The  thermocline  is 
centered  at  depth  z0  (=  -50  m)  with  thickness  D  (  =  120  m).  The 
buoyancy  frequency  derived  from  Eq.  (8)  is 

N  =  [-(g/p0)dpa/dz]l'2=  (1^5)  sech{^r)  =NmaxSec/j(?7p). 

0) 

The  maximum  buoyancy  frequency  Nmax  is  located  at  z=z0. 
For  A p  =  6kg /m3  and  p0  =  1028 kg/m3,  Nmax  =  N0  =  0.01 5 rad /s. 
Mode-1  waves  associated  with  Eq.  (9)  have  zero  horizontal  velocity 
and  maximum  vertical  velocity  at  z=  -200  m  (see  Fig.  1  in  Shaw 
et  al.,  2009).  The  phase  speed  of  the  mode-1  wave  is  1.52  m/s  for 
Nmax  =  N0,  and  is  proportional  to  Nmax  in  other  experiments.  For 
semidiurnal  tides,  the  internal  wave  is  evanescent  below  z=  -  700  m. 

The  baroclinic  velocity  iT  =  (u'.v'.w')  is  calculated  by  subtract¬ 
ing  the  depth-averaged  horizontal  velocity  from  the  total  velocity 
u.  The  perturbation  kinetic  energy  density  is  defined  by 

K  =  ip0(u/2  +  v/2+w/2), 


V  u  =  0. 


(10) 
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and  the  perturbation  potential  energy  density  is 
2p0N2' 


(11) 


The  internal  wave  energy  density  E  is  the  sum  of  Eqs.  (10) 
and  (11). 

The  internal  wave  energy  flux  is  p'u',  where  p'  is  the  perturbation 
pressure.  The  perturbation  pressure  p'  is  obtained  by  subtracting 
the  depth-averaged  hydrostatic  pressure  p  from  p  (Nash  et  al., 
2005): 

P'(*,z,t)  =  p(x,z,t)-p(x,t).  (12) 


The  depth  averaged  pressure  is  given  by 


P  = 


p\x,z\t)gdz' 


(13) 


where  zh  is  the  z-coordinate  of  the  bottom.  In  this  paper,  the  sum 
of  the  maximum  horizontal  energy  flux  on  the  two  sides  of  the 
ridge  at  x=  -xB  and  xB  represents  the  total  internal  wave  energy 
produced  by  the  ridge 

F=~  f  P'u'\x  =  XBdz+  [  P'u'\x=-xdz ■  (I4) 

JZb  Jzh 


Table  1 

List  of  experiments.  Each  experiment  is  carried  out  for  ridge  heights  /io  =  200,  400, 
and  600  m. 


Experiments 

I  (km) 

U0  (m/s) 

Nmax 

la 

15 

0.1 

N0 

lb 

30 

0.1 

No 

lc 

45 

0.1 

No 

Id 

60 

0.1 

N0 

2a 

30 

0.05 

No 

2b 

30 

0.2 

N0 

3a 

30 

0.1 

N0/2 

3b 

30 

0.1 

21/2N0 

4a 

15 

0.1 

No/2 

4b 

15 

0.1 

21/2N0 

3.  Results 

Three  sets  of  experiments  with  h0= 200,  400,  and  600  m  are 
carried  out  in  this  study.  Table  1  lists  the  topographic  width  L,  the 
amplitude  of  the  barotropic  tidal  velocity  U0,  and  the  maximum 
buoyancy  frequency  Nmax  for  each  h0. 

3A.  Internal  wave  generation 

Fig.  1  shows  the  distribution  of  the  baroclinic  velocity  a!  during 
the  ebb  phase  of  the  fourth  tidal  cycle  in  experiment  la  with 
h0 =400  m.  The  barotropic  tidal  current  is  zero  in  the  top  and 
bottom  panels  and  eastward  in  the  two  middle  panels.  In  Fig.  1, 
waves  emit  from  the  ridge  crest  and  split  into  two  nearly  symmetric 
paths  in  the  strongly  stratified  upper  ocean.  Over  the  ridge,  the 
baroclinic  velocity  uf  is  initially  negative  and  becomes  positive 
following  the  development  of  the  positive  barotropic  tidal  current. 
Thus,  a  wave  beam  of  positive  (negative)  u!  over  the  ridge  top 
originates  at  the  beginning  of  the  ebb  (flood)  tide  (Shaw  et  al.,  2009). 

Downward  wave  reflections  by  the  surface  occur  at  x=  ±26 
and  ±75  km.  At  x=  ±  50  km,  waves  are  reflected  upward  from 
z=-200m,  the  lower  boundary  of  the  thermocline.  Mode-1 
waves  appear  beyond  x=  +  100  km.  Intensification  of  the  wave 
front  with  negative  u'  leads  to  a  tidal  bore  at  x=-120km  in 
Fig.  la.  This  bore  can  be  traced  back  to  the  start  of  the  flood  tide  in 
the  second  tidal  cycle.  Tidal  bores  with  positive  u!  form  on  the 
east  side  of  the  ridge.  These  tidal  bores  eventually  develop  into 
rank-ordered  internal  solitary  waves  (Shaw  et  al.,  2009). 

Internal  wave  energy  density  at  z=  -5  m  for  experiment  la  is 
shown  on  the  x-t  plane  in  Fig.  2.  Maxima  are  at  x=  ±26  and 
±75  km,  locations  of  the  first  and  second  wave  reflections  from 
the  surface  in  Fig.  1.  Propagation  of  mode-1  hydrostatic  waves  is 
shown  as  bands  of  energy  maxima  on  each  side  of  the  ridge.  The 
stronger  ones  are  associated  with  negative  (positive)  u!  on  the 
left-  (right-)hand  side  of  the  ridge  in  Fig.  la.  The  speed  of 
propagation  is  close  to  the  mode-1  phase  velocity  of  1.52  m/s.  The 
energy  density  of  the  mode-1  wave  is  nearly  uniform  during 
propagation.  Small-scale  solitary  waves  are  shown  as  strings  of 
localized  maxima;  the  data  used  for  plotting  are  stored  hourly, 
insufficient  to  resolve  the  propagation  of  short  waves  (Shaw  et  al., 
2009). 
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Fig.  1.  Contour  plots  of  the  baroclinic  velocity  a!  (m/s)  in  experiment  la  with  1  =  15  km  and  /io=400  m  at  t=37.26  to  43.47  h  with  time  increments  of  2  h.  The  barotropic 
tidal  flow  is  to  the  right  in  the  two  middle  panels.  The  contour  interval  is  0.04  m/s. 
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Fig.  2.  Contour  plot  of  internal  wave  energy  density  (J/m3)  at  z=  -5  m  on  the  x-t  plane  for  experiment  la.  The  contour  interval  is  2  J/m3. 
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Fig.  3.  Internal  wave  energy  density  (J/m3)  for  I  equal  to  (a)  15,  (b)  30,  (c)  45,  and  (d)  60  km  (experiments  la-d)  with  ridge  height  at  400  m.  Averaging  is  performed  on 
hourly  data  over  three  tidal  cycles  between  t= 25.88  and  62.10  h.  The  contour  intervals  in  the  four  panels  are  2,  0.5,  0.2,  and  0.2  J/m3,  respectively. 


3.2.  Effects  of  topographic  width 

Fig.  3  compares  the  internal  wave  energy  density  for  ridges  of 
four  different  widths  (experiments  la-d)  with  the  same  ridge 
height  h0 =400  m.  The  energy  density  is  the  average  of  the  hourly 
data  over  3  tidal  cycles  from  25.875  to  62.1  h.  Fig.  3a  shows  the 
result  from  experiment  la  with  the  steepest  ridge  slope  (L=  15 
km).  Internal  wave  beams  similar  to  those  in  the  u'  plot  in  Fig.  1 
are  clearly  seen.  Between  200  and  400  m  depths,  wave  beams  are 
nearly  vertical  because  of  the  weak  stratification  in  the  water 
column.  In  the  strongly  stratified  upper  200  m  of  the  water 
column,  wave  beams  are  refracted  toward  the  horizontal  direc¬ 
tion.  Symmetric  wave  beams  appear  on  the  two  sides  of  the  ridge. 
Wave  reflections  from  the  surface  occur  near  the  edge  of 
topographic  influence  and  produce  energy  density  maxima  at 
the  surface.  After  reflection,  internal  waves  propagate  away  from 
the  ridge  in  an  upper  ocean  wave  guide. 


Increasing  L  to  30  km  reduces  the  energy  density  and  widens 
the  wave  beam  (Fig.  3b).  Internal  wave  beams  are  barely  visible, 
but  maxima  in  energy  density  are  still  at  the  surface.  At  1=45  and 
60  km,  distribution  of  the  energy  density  differs  greatly  from  that 
for  a  narrow  ridge  (Fig.  3c, d).  The  maximum  energy  density  is 
much  smaller  than  that  in  Fig.  3a.  Furthermore,  reflections  from 
the  surface  produce  maxima  at  z=  - 125  m  instead  of  the  surface. 
Fig.  3  suggests  that  the  ridge  with  L= 15  km  is  supercritical,  and 
the  ridges  with  1=45  and  60  km  are  subcritical.  The  slope  is  near 
the  critical  value  for  the  ridge  with  1=30  km.  A  supercritical  ridge 
is  most  efficient  in  producing  internal  waves. 

The  vertical  energy  flux  w'p'  at  t=41.40  h  is  shown  in  Fig.  4  for 
experiments  la  and  Id  with  1=15  and  60  km,  respectively.  Over  a 
narrow  ridge  (Fig.  4a),  energy  flux  between  z=  -  700  and  -  200  m 
is  large  and  upward  on  the  flank  of  the  ridge.  Above  z=  -200  m, 
upward  and  downward  fluxes  follow  the  path  of  the  wave  beam 
shown  in  Fig.  3a.  Free  propagating  hydrostatic  waves  occur  at 
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Fig.  4.  Distribution  of  the  vertical  energy  flux  (W/m2)  at  41.40  h  in  experiments  (a)  la  and  (b)  Id  with  ridge  height  at  400  m.  The  contour  intervals  are  0.02  and  0.002 
W/m2,  respectively. 


t  =  25.875-62.1  hours 


Fig.  5.  Vertically  integrated  horizontal  energy  flux  (W/m)  as  a  function  of  x  for  L=  15, 30, 45,  and  60  km  with  ridge  height  at  400  m  (experiments  la-d).  The  flux  is  averaged 
over  three  tidal  cycles  from  25.88  to  62.1  h.  The  ridge  is  located  at  x=0  km. 


|x|  >60  km.  When  L= 60  km  (Fig.  4b),  energy  flux  is  mainly  over 
the  flank  of  the  ridge,  upward  below  z=  -200  m  and  downward 
in  the  upper  200  m  similar  to  a  mode-1  standing  wave.  The  flux  is 
nearly  10  times  smaller  than  in  experiment  la,  and  little  energy 
leaves  the  ridge  topography. 

Fig.  5  shows  the  distribution  of  vertically  integrated  horizontal 
energy  flux  u'p'dz  as  a  function  of  x  in  experiments  la- Id.  The 
flux  is  averaged  over  3  tidal  periods  as  in  Fig.  3.  The  horizontal 
energy  fluxes  increase  rapidly  with  x  near  the  ridge  crest,  reach 
maximum  values  at  the  location  of  maximum  energy  density  at 
the  surface,  and  decay  slightly  thereafter  (Fig.  3).  The  maximum 
horizontal  energy  flux  is  highest  in  experiments  using  the 
narrowest  ridge  (L= 15  km)  and  decreases  rapidly  with  increasing 
ridge  width.  The  ridge  half-width  gives  the  proper  horizontal 
length  scale  in  this  plot. 

3.3.  Effects  of  the  barotropic  tidal  current 

Fig.  6  compares  the  internal  wave  energy  density  generated  by 
barotropic  tidal  currents  of  amplitudes  U0= 0.2,  0.1,  and  0.05  m/s 
over  a  400-m  ridge  with  L= 30  km.  The  pattern  of  internal  wave 
energy  density  is  qualitatively  similar  in  the  three  experiments; 
only  the  magnitudes  are  different.  Ripples  in  the  experiment  of 
the  strongest  tidal  current  indicate  the  formation  of  internal 
solitary  waves  (Fig.  2). 

Fig.  7  shows  the  vertically  integrated  horizontal  energy  flux 
normalized  by  Uo  for  the  three  experiments  plotted  in  Fig.  6.  After 
normalization,  the  three  curves  become  similar,  especially  in  the 
wave  generation  region.  Fig.  7  confirms  that  the  internal  wave 


energy  flux  is  proportional  to  the  square  of  the  amplitude  of  the 
barotropic  tidal  current  as  suggested  in  the  linear  theory.  In  the 
range  of  parameters  used  in  this  study,  varying  the  strength  of  the 
barotropic  tidal  current  does  not  affect  the  generation  of  internal 
waves  qualitatively. 


3.4.  Effects  of  stratification 

The  perturbation  density  fields  at  37.26  h  in  experiments  3b  and 
3a  with  Nmax=21/2N0  and  N0/2,  respectively,  are  shown  in  Fig.  8.  The 
ridge  height  h0  is  400  m  and  L= 30  km.  When  Nmax=21/2N0,  wave 
beams  form  over  the  ridge  crest  (Fig.  8a).  The  angle  of  the  wave 
beam  with  the  vertical  is  greater  than  that  at  Nmax=N0  in  Fig.  4a, 
consistent  with  the  internal  wave  theory.  Large-amplitude  internal 
waves  of  wavelength  greater  than  the  scale  of  topography  propagate 
away  from  the  ridge.  In  weak  stratification  (Fig.  8b),  internal  waves 
generated  over  the  ridge  are  weak  and  form  a  mode-1  standing 
wave  over  the  ridge  topography.  Patterns  of  the  vertical  energy  flux 
in  experiments  3b  and  3a  are  similar  to  those  in  Fig.  4  (not  shown). 

Fig.  9  shows  the  average  energy  density  of  internal  waves  in 
experiments  3b,  lb,  and  3a  with  maximum  buoyancy  frequency 
Nmax=21/2N0,  N0,  and  N0/2,  respectively.  The  ridge  height  h0  is 
400  m  and  L= 30  km.  The  distribution  of  energy  density  differs 
qualitatively  in  these  experiments.  As  shown  in  Fig.  3b,  the  ridge 
slope  in  experiment  lb  is  near  the  critical  value  (Fig.  9b).  When 
Nmax=2 1/2N0,  vertically  propagating  wave  beams  in  the  upper 
ocean  are  distinct  (Fig.  9a).  On  the  other  hand,  standing  waves 
form  over  the  ridge  for  Nmax=N0/2  (Fig.  9c),  similar  to  the 
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Fig.  6.  Distribution  of  internal  wave  energy  density  (J/m3)  for  U0  equal  to  (a)  0.2,  (b)  0.1,  and  (c)  0.05  m/s  with  h0=400  m  and  L= 30  km  (experiments  2b,  lb,  and  2a).  The 
energy  density  is  averaged  from  25.88  to  62.1  h.  The  contour  intervals  are  2,  0.5,  and  0.1  J/m3,  respectively. 


x  (km) 

Fig.  7.  Vertically  integrated  horizontal  energy  flux  (W/m)  as  a  function  of  x  for  the  three  experiments  shown  in  Fig.  6.  The  energy  flux  is  normalized  by  Do  and  averaged 
over  three  tidal  cycles  from  25.88  to  62.1  h. 


experiment  with  L=60km  in  Fig.  3d.  The  maxima  in  energy 
density  in  Fig.  9c  are  located  below  surface.  The  energy  flux  of  the 
internal  wave  is  largest  in  Fig.  9a.  Increasing  stratification  has  the 


same  effect  on  internal  wave  generation  as  increasing  the  ridge 
slope.  When  stratification  is  strong,  wave  beams  form  and  a  large 
amount  of  internal  wave  energy  radiates  away  from  the  ridge. 
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Fig.  9.  Distribution  of  internal  wave  energy  density  (J/m3)  for  ridges  at  ho=400  m  and  L  =  30  km  in  experiments  (a)  3b,  (b)  lb  and  (c)  3a.  Time  average  is  over  the  same 
period  as  in  Fig.  3.  The  contour  intervals  are  1,  0.5,  and  0.1  J/m3,  respectively. 
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Fig.  10.  Distribution  of  internal  wave  energy  density  (J/m3)  for  varying  ridge  heights  at  L=  15  km.  The  energy  density  is  averaged  over  three  tidal  cycles  from  25.88  to 
62.1  h.  The  contour  intervals  are  10,  2,  and  0.2  J/m3,  respectively. 


3.5.  Effects  of  ridge  height 

Fig.  10  shows  the  energy  density  of  the  internal  wave  for  a 
steep  ridge  with  1=15  km  at  h0= 200,  400,  and  600  m.  The 
parameters  are  the  same  as  those  in  experiment  la.  The  wave 
patterns  in  these  experiments  are  similar.  Above  z=-200m, 
wave  beams  form  at  nearly  the  same  locations  in  all  three 
experiments.  However,  the  maximum  energy  is  much  reduced 
when  the  ridge  height  decreases.  The  ratio  of  energy  density  from 
the  smallest  ridge  to  the  tallest  ridge  is  1:7.5:30. 

Varying  the  ridge  height  h0  in  experiment  Id  (L= 60  km)  does 
not  change  the  pattern  of  the  internal  wave;  subsurface  maxima 
appear  over  the  ridge  in  all  experiments  (Fig.  11).  Similar  to  the 
experiments  with  a  narrow  ridge  in  Fig.  10,  maximum  energy 
density  decreases  rapidly  with  decreasing  ridge  height.  The  energy 
ratio  for  the  three  ridges  is  1:2:20.  In  the  linear  theory  (Llewellyn 
Smith  and  Young,  2002),  the  ratio  is  1:5:24,  the  buoyancy 
frequency  ratio  at  the  depth  of  the  ridge  crest.  The  energy  ratios 


in  the  numerical  experiments  have  the  same  order  of  magnitude  as 
that  obtained  from  the  linear  theory  but  the  linear  theory 
underestimates  (overestimates)  the  energy  density  of  the  internal 
wave  when  the  topographic  slope  is  supercritical  (subcritical). 


4.  Parameterization  of  the  energy  flux 

Experiments  using  parameters  shown  in  Table  1  are  repeated 
for  three  ridge  heights  h0= 200,  400,  and  600  m.  From  these 
experiments,  a  parameterization  scheme  can  be  formulated  to 
describe  the  dependence  of  the  vertically  integrated  horizontal 
energy  flux  on  the  ridge  slope,  barotropic  tidal  current  amplitude, 
and  stratification  for  each  topographic  height.  Normalization  of 
the  energy  flux  based  on  the  linear  theory  is  modified  to  fit  the 
numerical  data  for  tall  ridges.  The  scheme  is  represented  by  a 
nondimensional  relationship  between  the  normalized  energy  flux 
and  the  slope  parameter  y. 
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Fig.  11.  Same  as  Fig.  10  except  for  1  =  60  km.  The  contour  intervals  are  2,  0.1,  and  0.05  J/m3,  respectively. 


4A.  Normalization  scheme 


The  experiments  in  Section  3.3  show  that  U$  is  an  appropriate 
factor  for  energy  flux  generated  over  finite  ridges.  Physically,  the 
energy  flux  should  be  proportional  to  the  depth-integrated 
barotropic  tidal  energy.  Thus,  it  is  reasonable  to  scale  the  energy 
flux  by  HU%.  In  the  linear  theory,  the  buoyancy  frequency  at  the 
ridge  crest  provides  another  scaling  factor.  However,  Fig.  4  shows 
that  the  vertical  energy  flux  is  affected  by  a  large  area  of  the  ridge 
flank,  not  just  the  ridge  crest.  A  better  scale  would  be  the  averaged 
buoyancy  frequency  over  the  depths  of  the  ridge 


/ 


-H  +  ft0 
-H 


Ndz. 


(15) 


Using  Eq.  (15),  the  local  Froude  number  becomes 


Nrhom 


(16) 


The  normalized  energy  flux  is  defined  as 


pHU02Nrh0 ' 


(17) 


If  multiplied  by  <5,  the  normalization  factor  in  Eq.  (17)  is  the  same 
as  that  of  Eq.  (1)  for  N^>co  and  /=  0  in  an  ocean  of  constant 
stratification. 

The  slope  parameter  y  depends  on  the  ridge  slope  and  the 
slope  of  the  wave  beam;  both  are  not  constant  in  this  study.  In 
Section  3.5,  raising  or  lowering  the  ridge  affects  the  energy 
density  of  the  internal  wave  but  not  the  wave  beam  in  the  upper 
200  m  of  the  water  column  (Figs.  10  and  11).  Therefore,  the 
maximum  ridge  slope,  which  does  not  vary  in  these  experiments, 
is  a  good  description  of  y.  Furthermore,  Fig.  4a  shows  that  tilting 
of  the  wave  beam  in  the  strongly  stratified  upper-ocean  wave 
guide  could  be  represented  by  the  averaged  wave  slope  between 
zt=-  200  m  and  0  m 


s 


1/2 

dz. 


(18) 


Table  2 

The  nondimensional  horizontal  energy  flux  (F),  the  slope  parameter  (y),  and  the 
local  Froude  number  (Fr)  calculated  from  all  experiments  in  this  study. 


Ridge  height  (h0) 

experiment 

600  m 

400  m 

200  m 

F' 

y 

Fr 

F' 

y 

Fr 

F' 

y 

Fr 

la 

1.16 

2.6 

0.1 

1.35 

2.6 

0.56 

0.62 

2.6 

2.5 

lb 

0.81 

1.3 

0.1 

1.01 

1.3 

0.56 

0.81 

1.3 

2.5 

lc 

0.30 

0.87 

0.1 

0.33 

0.87 

0.56 

0.74 

0.87 

2.5 

Id 

0.06 

0.65 

0.1 

0.12 

0.65 

0.56 

0.35 

0.65 

2.5 

2a 

0.84 

1.3 

0.05 

1.02 

1.3 

0.28 

0.81 

1.3 

1.25 

2b 

0.74 

1.3 

0.2 

0.96 

1.3 

1.11 

0.80 

1.3 

5 

3a 

0.05 

0.65 

0.2 

0.12 

0.65 

1.19 

0.32 

0.65 

5 

3b 

1.08 

1.83 

0.07 

1.24 

1.83 

0.37 

0.80 

1.83 

2 

4a 

0.70 

1.3 

0.2 

1.07 

1.3 

1.19 

0.83 

1.3 

5 

4b 

1.40 

3.67 

0.07 

1.12 

3.67 

0.37 

0.58 

3.67 

2 

The  slope  parameter  y  is  the  ratio  of  the  maximum  topographic 
slope  and  s. 

Table  2  summarizes  the  nondimensional  parameters  obtained 
from  all  experiments  in  this  study.  The  Froude  number  is  less  than 
1.2  for  h0 =400  and  600  m  and  is  greater  than  1.2  for  h0= 200  m. 
The  normalized  horizontal  energy  flux  F  is  plotted  in  Fig.  12  for 
these  two  groups.  Two  curves  are  drawn  to  describe  the 
dependence  of  F  on  y.  In  Fig.  12a,  F  is  nearly  constant  for  large 
y  and  decreases  to  zero  rapidly  at  y<  1.  The  transition  occurs 
between  y  =  0.5  and  1.5.  It  seems  that  energy  conversion  by  wave 
beams  is  possible  only  when  the  ridge  slope  is  supercritical. 
Similarly,  the  normalized  energy  flux  in  Fig.  12b  is  nearly  zero  for 
subcritical  ridges  as  in  Fig.  12a.  In  the  supercritical  range,  F 
decreases  slightly  with  increasing  y.  Fig.  12  suggests  that 
barotropic-to-baroclinic  energy  conversion  occurs  only  when  y 
is  greater  than  1. 

4.2.  Comparison  with  observations 

Table  3  lists  the  energy  fluxes  of  strong  internal  waves  at  four 
locations  in  the  world  ocean  from  earlier  observations  and 
numerical  studies.  The  energy  flux  in  the  Luzon  Strait  is  from 
the  real-time  simulation  at  Naval  Research  Laboratory  for  2005 
(S.  F.  Lin,  2009,  personal  communication).  At  French  Frigate,  the 
maximum  internal  wave  energy  flux  is  39  kW/m  from  velocity 
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a 
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Fig.  12.  Plots  of  normalized  horizontal  energy  flux  (F)  in  Table  2  as  a  function  of 
the  slope  parameter  (y):  (a)  for  ridge  heights  at  400  m  (triangles)  and  600  m 
(circles)  and  (b)  for  ridge  heights  at  200  m  (crosses).  In  the  top  panel,  *  represents 
estimates  in  the  ocean  at  Bay  of  Biscay  (B),  Luzon  Strait  (L),  French  Frigate  (F)  and 
Kauai  Channel  (I<)  from  Table  3.  The  three  regimes  for  small,  near  critical,  and  large 
y  are  enclosed  by  a  rectangle,  triangle,  and  ellipse,  respectively. 


Table  3 

Energy  fluxes  (F),  depth  range  of  ridges,  the  slope  parameter  (y),  and  the  amplitude 
of  barotropic  tidal  current  ( U0 )  at  locations  where  strong  internal  waves  are 
observed. 


Location 

F  (kW/m) 

Depth  range  (m) 

7 

U0  (m/s) 

Luzon  Strait 

55 

500-3000 

3.42 

0.06 

French  Frigate 

40 

900-3400 

3.93 

0.05 

Kauai  Channel 

26 

1000-3000 

3.61 

0.05 

Bay  of  Biscay 

10 

250-4000 

2.14 

0.025 

measurements  (Rudnick  et  al.,  2003)  and  45  kW/m  from  internal 
tide  simulation  (Simmons  et  al.,  2004).  The  averaged  value  is 
included  in  Table  3.  The  energy  flux  in  Kauai  Channel  is  from  the 
estimate  of  Simmons  et  al.  (2004).  The  energy  flux  at  Bay  of  Biscay 
is  from  the  current  meter  measurements  of  Gerkema  et  al.  (2004). 

The  bottom  depth  is  obtained  from  the  Gridded  Global  Relief 
Data  (ETOP02)  with  2-min  resolution  (http://www.ngdc.noaa. 
gov/mgg/fliers/01  mgg04.html).  Topography  is  smoothed  by  a  2-D 
Savitzky-Golay  filter  (Ratzlaff  and  Johnson,  1989;  Kuo  et  al., 
1991).  The  topographic  slope  is  then  calculated  from  the 
smoothed  topography  as  the  first  derivatives  with  respect  to  x 
and  y.  The  buoyancy  frequency  is  calculated  from  the  climatolo¬ 
gical  temperature  and  salinity  fields  of  Levitus  (1998,  http://www. 
esrl.noaa.gov/psd/data/gridded/data.nodc.woa98.html).  The  ampli¬ 
tude  of  the  barotropic  tidal  current  is  obtained  from  previous 
studies  (e.g.  Legg  and  Klymak,  2008;  Gerkema  et  al.,  2004;  Jan  et  al., 
2008). 

The  depth  range  of  the  ridge  is  obtained  from  the  smoothed 
topography.  The  ridge  slope  is  the  average  topographic  slope  in 
this  range.  The  wave  slope  is  estimated  using  Eq.  (18).  Large 
uncertainties  exist  in  the  calculation  of  maximum  topographic 
slope  and  thusy.  Nevertheless,  the  estimated  y  is  greater  than  1  at 
all  four  locations.  The  vertical  integral  of  the  buoyancy  frequency 
in  the  depth  range  of  the  ridge  is  obtained  from  Eq.  (15)  for 
normalization  of  the  energy  flux  and  calculation  of  the  local 
Froude  number  (16).  The  Froude  number  is  much  smaller  than  1 
at  all  locations.  The  normalized  energy  fluxes  are  nearly  constant 
for  the  four  locations  and  follow  the  curve  in  Fig.  12a. 


5.  Discussion 

In  this  study,  the  wave  beam  is  nearly  vertical  below  200  m 
because  of  weak  stratification  in  the  deep  water.  The  slope 
parameter  y  for  the  supercritical  ridges  at  200  and  400  m  would 
be  much  smaller  than  1  if  calculated  from  the  depth-averaged 
buoyancy  frequency  as  suggested  by  Llewellyn  Smith  and  Young 
(2002).  Figs.  10  and  11  validate  the  use  of  maximum  wave  slope 
in  the  calculation  of  y.  Consequently,  a  ridge  could  become 
supercritical  at  the  presence  of  a  strong  upper-ocean  thermocline 
even  when  the  stratification  at  the  ridge  crest  is  weak.  The  critical 
y  obtained  in  this  study  is  about  1,  independent  of  the  ridge 
height.  The  regime  of  strong  wave  beams  with  y  >  1.5  is  grouped 
in  an  ellipse  in  Fig.  12.  The  wave  regime  with  y  <  1  is  enclosed  by 
a  rectangle.  In  the  latter  experiments,  no  internal  wave  beams 
form,  and  internal  waves  of  vertical  modes  develop  over  the  ridge. 
The  experiments  enclosed  by  triangles  show  transition  from  a 
modal  structure  to  wave  beams. 

Upward  baroclinic  energy  flux  is  produced  on  the  flank  of  the 
ridge  by  the  tidal  current.  On  a  wide  ridge  (y  <  1 ),  downward 
propagating  waves  from  reflection  of  the  surface  interfere  with 
the  upward  propagating  waves  over  the  ridge,  producing  a  mode- 
1  standing  wave.  The  vertical  mode  does  not  match  that  of  free 
propagating  waves  in  the  undisturbed  ocean;  little  energy 
propagates  away  from  the  ridge.  Over  a  supercritical  ridge 
(y  >  1 ),  surface  reflection  of  the  wave  beam  occurs  outside  ridge 
topography.  Not  influenced  by  the  upward  flux  on  the  flank  of  the 
ridge,  narrow  slanting  wave  beams  form  in  the  strongly  stratified 
upper  ocean  between  the  surface  and  the  base  of  the  thermocline. 
The  wave  beams  then  develop  into  free  propagating  hydrostatic 
waves  with  vertical  structure  matching  that  in  the  undisturbed 
ocean.  A  large  amount  of  baroclinic  energy  radiates  away  from  the 
ridge. 

The  magnitude  of  the  energy  flux  is  proportional  to  HUq  and 
the  integral  of  N  over  the  depths  of  the  ridge.  Over  a  supercritical 
ridge,  the  normalized  energy  reaches  a  constant  value  for  ridge 
heights  at  400  and  600  m  (Fig.  12a).  The  local  Froude  number  is 
small  in  these  experiments  (Table  2).  When  Fr<l  and  y>l, 
overturning  may  limit  the  growth  of  isopycnal  displacement  (Legg 
and  Klymak,  2008),  resulting  in  saturation  of  the  energy  flux  for 
large  y.  The  decrease  in  energy  flux  with  y  in  the  supercritical 
range  for  200-m  ridges  in  Fig.  12b  is  likely  an  artifact  of  lowering 
the  ridge  crest  close  to  bottom.  The  small  surface  area  of  the 
topography  reduces  the  total  vertical  flux  generated  on  the  flank 
of  the  ridge. 

In  an  ocean  with  depth-varying  buoyancy  frequency,  the  WKB 
approximation  gives  the  scaled  ridge  height 

B  nJZ^Ndz  nhpNr 
f°H  Ndz  HN 

where  N  =  ^  f°H  Ndz  and  Nr  given  by  (1 5).  In  terms  of  B,  Llewellyn 
Smith  and  Young  (2003)  provided  a  formula  for  energy  conversion 
from  a  knife-edge  ridge.  In  the  case/=0,  the  flux  is 

Fk„ife  =  p0U20NrH2  =  PoHU^Nrho^^LMiB/n), 

n  4HN 

where  M(B/n)  =  (4I<(B/n)/nB2).  The  function  M  for  density 
stratification  (9)  is  between  1  and  1.03  if  the  ridge  height  is 
below  600  m.  The  buoyancy  frequency  at  the  ridge  crest  is  shown 
in  Fig.  13a.  After  normalization  by  Eq.  (17),  the  energy  flux  for  a 
knife-edge  ridge  is  F'knife  =  (57iNrNr/4N  ),  which  is  plotted  in 
Fig.  13b.  The  flux  at  fro=600m  is  about  half  that  over  a 
supercritical  bell-shaped  ridge  of  the  same  height  (Fig.  12a). 
While  the  normalized  flux  in  Fig.  12a  is  independent  of  ridge 
height,  F  calculated  from  the  linear  theory  for  a  knife-edge  ridge 
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Fig.  13.  (a)  Buoyancy  frequency  at  the  ridge  crest  (Nr)  from  (9)  as  a  function  of  the 
ridge  height  (h0)  for  a  knife-edge  ridge  with  No=0.015  rad/s.  (b)  The  normalized 
conversion  rate  (F)  calculated  using  the  formula  of  Llewellyn  Smith  and  Young 
(2003)  after  normalization  by  (17). 


is  very  sensitive  to  the  ridge  height.  In  the  analytical  solution,  F 
for  a  400-m  ridge  is  about  4%  of  that  at  600  m.  Since  wave  beam 
formation  in  the  upper  thermocline  governs  energy  conversion  in 
this  study,  it  is  not  surprising  that  the  normalized  flux  in  Fig.  12a 
is  independent  of  stratification  at  the  ridge  crest. 


When  the  Froude  number  is  small,  a  supercritical  ridge  produces 
maximum  energy  flux  independent  of  the  slope  parameter.  Since 
the  energy  flux  is  proportional  to  the  stratification  at  depths  of  the 
ridge  according  to  Eq.  (17),  only  tall  ridges  reaching  strongly 
stratified  water  contribute  significantly  to  the  energy  flux  in  the 
world  ocean.  Furthermore,  the  normalized  energy  flux  produced 
by  a  tall  supercritical  ridge  is  insensitive  to  ridge  height  and  the 
ridge  slope.  Fig.  12a  suggests  a  practical  method  to  estimate  the 
energy  flux  of  internal  waves  in  the  world  ocean. 

In  the  northern  South  China  Sea,  energy  of  the  internal  wave 
increases  during  the  spring  tide  according  to  Eq.  (17).  In  addition, 
energy  flux  may  become  large  when  the  Kuroshio  front  increases 
the  stratification  at  the  ridge  crest  in  the  Luzon  Strait.  More 
importantly,  the  slope  parameter  may  become  supercritical  at  the 
presence  of  a  strong,  shallow  thermocline  in  the  Luzon  Strait. 
The  development  of  a  strong  seasonal  thermocline  may  explain 
the  frequent  occurrences  of  internal  solitary  waves  in  the 
northern  South  China  Sea  in  spring  and  summer. 
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